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II. 1 GEOMETRY 


The gross geometric characteristics of the aircraft under study are 
specified to the subroutine SIZE which follows the sequence shown in Figure 
II. 1.1. The principal quantities specified are both geometric (lengths and 
areas) and operational (altitude, Mach number). The sequence of computations 
carried out by SIZE is controlled by the parameter NPC which is passed into 
SIZE by COMMON along with the input quantities. When NPC«0 the computation 
is initialized. Subsequently, NPC is set to 2, and the program advances 
through the geometric computations following the sequence shown in Figure 
II. 1.1. 

II. 1.1 Fuselage Geometry 

The fuselage is assumed to consist of four components as illustrated in 
Figure II. 1.2: a nose cone, a pilot's compartment, a cylindrical cabin, 

and a tail cone. The external cabin width is found in inches as 


* W + A * 

C AB S S 


W 


AS 


+ 12 


where 

S * number of passenger seats abreast 

AB 

W g ■ seat widths, in 

A ■ number jf aisles 

W, - aisle width, in 

AS 


(II. 1.1) 


and where the cabin walls are assumed to six inches thick on each side. This 
dimension is next converted to feet, 


H 


C 


S 


WF 


V 12 


(II. I..**) 
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FIGURE II. 1.1 


SUBROUTINE SIZE, COMPUTATIONAL FLOW 




and the cabin length is estimated in feet as 


L C ■ <P «C - 1 ,P S /(12 V 


(II. 1.3) 


where 


■ number of passengers - x 
P “ seat pitch, in. 


Next, the nose cone diameter is found as the difference of the body 
diameter and the windshield height. 


«N " H C - H CK 


(II. 1.4) 


as shown in Figure II. 1.2, and if S AB — 2, the cabin computations are complete. 
If S ( 2, (in-line seating) however, the cabin length is found by 

AB 


L c “ P AX * V 12 ' 


(II. 1.5) 


and the nose and tail cone diameters are 


h n - s 


WF 


H - H + H 
C N CK 


(11. 1.6) 

(11. 1.7) 


where H is the input windshield height. 

The total fuselage length is then given by the sum 


EL - ELOD * H + EL__ + L_ + ELOD_ * H 
F N N PC C T C 


(II. 1.8) 


where 


ELOD^ - length to diameter of nose cone 
EL pc • length of pilot's cockpit 
ELOD • length to diameter of tail cone 
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The nttmil fuselage watted area ia next given by 


S P " H c I 2 * 5 < ELOD N * Hj, ♦ “-pc) ♦ 3.14 L c ♦ 2.1 BLOO t * H C 1 

(II. 1.9) 

where the first and last numerical coefficients approximate the effects of the 
nose and tail taper ratios. If the configuration incorporates a tail boom 
instead of a cone ( K C0 NFG B ^ ' t * ien t * ie boom length is 


where 


“*M " EL F - ’SfC 

^FPC “ ^ D N * »N + EL PC + L C 


(II. 1.10) 


(II. 1.11) 


The fuselage wetted area in this case does not include the tail cone portion, 
i . e . , 

S P * H c [2.5(ELOD n * H n + EL^.) + 3.14 * 1^] (II. 1.12) 

and the boom area is found assuming it to be cylindrical in shape: 


S WB* ■ 3 - 14 V%D 


(II. 1.13) 


where " tail boom length to diameter ratio. The cross sectional area 

of the tail boom is then given as 


“V ■ W“- * “Yod 


(II. 1.14) 


II. 1.2 Wing Geometry 
The wing area and span are given by 


s w - V"« 


(II. 1.15) 
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B - **•*„ 

where 

W Q ■ aircraft gross weight, Id. 
w ■ wing loading of aircraft, lb per aq ft 

CiS 

AR ■ aspect ratio of wing 

Wing root chord at the fuselage .enter line is obtained frcn 


(II. 1.16) 


CR CUI ■ 2 V lB,1 *V cu.1.17) 

and the mean aerodynamic chord, 

CBAR W • (2/3)CR CLW ll ♦ SuI - S^U * s^l (II. 1.18) 

where 

S - taper ratio of wing (tip chord/root chord) 

UN 

The leading edge sweep angle is found as a function of the input quarter 
chord sweep angle (DLM^) 

T SWPLE • t ' m<Sm V • (l - S U, >/[Wi(1 ♦ V’ * t "’ lDLM C4 ) 

(II. 1.19) 

and a similar expression gives the trailing edge sweep as 


T SWPTt ■ “"‘ S “V " 3 (S U) - ll/l ‘* <1 * V * “"‘““W 

(II. 1.20) 


The exposed root chord of the wing is next found as 


- CR - ♦ r tm * T / 2 

ROOTW CLW P KP SVFEPTE 


(XI. 1.21) 


where 


H 


P 


HP 


T /2 

SWPUC 


(II. 1.22) 
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V 


(XX. 1.23) 


as shown In Figure XI. 1.3. 

II. 1.3 Empennage Geometry 

Preliminary tail sizing is performed in SIZE based on tail volume coeffi- 
cients. Subsequent tail sizing may be performed in subroutine WGHT and/or 
TAIL according to the tail sizing option (LCWING) . 

Volume coefficients for the vertical and horizontal tails may be input as 
VBARVX and VBARHX respectively. Otherwise they are determined as empirical 
functions of fuselage and wing geometry as illustrated in Figure II. 1.4: 

V*H ■ ■ t ’ 5 E S S WF 2/(S W C “V * Si 



S * 0, for low tail 
AH 

S » 1, for T tail 
AH 


The horizontal and vertical tail areas are then computed from the volume 
coefficients: 


S - V SC 
HT BARH W H2 


(II. 1.27) 
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(II. 1.28) 


where 


S - V SC 
VT BAKV W V2 


WING MAC _ 71 

C H2 “ horizontal tail moment aim “ * 271 * 0955 S AH 


(II. 1.29) 


or input as COELTH 


°V2 " 


wing span 


vertical tail moment aim 


- 1.862 + .338 S 


AH 


(II. 1.30) 


or input as BOELTV 


The horizontal tail area is constrained to be less than or equal to S^/2; 
if as computed above exceeds this limit, both tail volume coefficients 
are re-computed : 


V BARH * ‘ 5 ^ C H2 


(II. 1.31) 


V ■ 1 V 

BAKV * BARH 


(II. 1.72) 


New tail areas are then computed based on these volume coefficients. 
The horizontal and vertical tail moment arms are computed from 


EL TH " V BARH S W C BARW^ S HT 


EL - V S B/S 

TV BAKV W VT 


Once the tail areas are determined, various geometric characteristics of each 
tail surface are computed in the same manner as was done for the wing. For 
example, for the horizontal tail 

1/2 

Span! B HT " (S HT *V (II. 1.33) 
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root chord t 


(II. 1.34) 


“cLHT ' 2 * 

■ran ««ro chord* CR CLHT * 2 * S ht /1B ht I1+SL Mh >1 (II. 1.35) 

where SL^ Is the horizontal tail taper ratio. Identical equations apply to 
the vertical tail. 


11.1.4 Cabin Pressurization 

The subroutine TPALT is called twice to give static pressures inside and 

outside the fuselage cabin. The inside pressure is specified by the cabin 

altitude of 8000 feet and the outside by the cruise altitude, HNCRU , and these 

altitudes are associated respectively with the pressures _ and P. ... 

CAB AMBC 

The input pressure differential DEL p is then checked to see if it is suffi- 
cient to maintain the 8000 ft cabin at the design cruise altitude; if not, 
the proper DEL p is then calculated. 

I I. 1.5 Nacelle Area 

Finally, if the engine nacelle size is input, " 2, the total 

nacelle wetted area is computed in square feet as 


S - EN_ * PI 
N P 


• D 


* 

BARN 


EL 

N 


where 

EN p ■ number of engines 

D - average nacelle diameter, ft 

BARN 

EL^ ■ nacelle length, ft 

PI - 3.1416 


(II. 1.36) 
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II. 1.6 Geometry Summary 


The calculations carried out by subroutine SIZE have defined the vehicle 
geometry with sufficient detail to permit the design to proceed to definition 
of aerodynamics, propulsive, weights, performance and economic characteristics. 
Section II. 2 of this report defines all input parameters and output charac- 
teristics of subroutine SIZE together with the computer program symbology 
employed . 
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II. 2 GEOMETRY MODEL USER’S MANUAL 

Subroutine SI2E is a relatively short program (110 cards) which requires 
about 35 input parameters related to the aircraft geometry and the design 
flight condition. These input variables are passed to SIZE from the MAIN 
control program through COMMON. The 36 output parameters are principally 
lengths and areas calculated according to the equations discussed in earlier 
sections of this volume. Most of the input data are specified by the user, 
since this is the first subroutine called by program MAIN, an d no other sub- 
routines have yet generated numerical results. However, many of the output 
quantities of SIZE will be needed as input to subsequent subroutines. 

In addition to the input parameters two internally generated control 
parameters are employed in SIZE. These are NPC which determines the SIZE 
function as follows: 

NPC - 1, 2 Geometry Computation 

■ 3 Geometry Initialization 

and ALTCAB * Cabin design altitude (Nominally set 

to 8000 feet in data statement 

A list of input parameters to SIZE is presented in Figure II. 2.1. Character- 
istics computed by subroutine SIZE are tabulated in Figure II. 2. 2. 
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FIGURE II. 2.1 SUBROUTINE SIZE INPUT PARAMETERS 


VARIABLE 

DESCRIPTION 

AR 

aspect ratio of wing 

ARHT 

aspect ratio of horizontal tail 

ARVT 

aspect ratio of vertical tail 

AS 

number of aisles 

BMLOD 

ratio of boom length to diameter 

BOELTV 

ratio of wing span to vertical tail moment arm 

COELTH 

ratio of wing chord to horizontal tail moment arm 

DBARN 

average nacelle diameter, ft 

DELP 

pressure difference, cabin minus ambient, lb per sq ft 

DLMC4 

quarter chord sweep angle, deg. 

ELN 

nacelle length, ft 

ELODN 

ratio of length to diameter of nose cone ' 

ELODT 

ratio of length to diameter of tail cone 

ELPC 

length of pilot* s cockpit 

ENP 

number of engine phylons 

HCK 

cockpit height, ft 

HNCRU 

cruise altitude, ft 

KCONFG 

0, standard fuselage 

1, tail boom fuselage 

KNAC 

0, nacelle drag computed as penalty to engine performance 
(turbofan only 

1, nacelle drag part of aerodynamic drag; nacelle sized 
by engine 

2, same as 1 except nacelle size input as DBARN, ELN 

PAX 

number of passengers 


1 1 - 2 
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FIGURE II. 2.1 SUBROUTINE SIZE INPUT PARAMETERS (Continued) 


VARIABLE 

DESCRIPTION 




PS seat pitch, inches 

SAB number of passenger seats abreast 

SAH dimensionless measure of horizontal tail height 

SLM taper ratio of wing 

SLMH taper ratio of horizontal tail 

SLMV taper ratio of vertical tail 

TCR, TCT thickness to chord ratios at root and tip of wing 

VBARHX horizontal tail sizing parameter 

VBARVX vertical tail sizing parameter 

WAS aisle width, inches 

WG gross weight, lb 

WGS wing loading, lb per sq ft 

WS seat width, inches 


1 1-2 
3 






FIGURE II. 2. 2 SUBROUTINE SIZE - OUTPUT CHARACTERISTICS 
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FIGURE II. 2. 2 SUBROUTINE SIZE - OUTPUT CHARACTERISTICS 
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II -3 GEOMETRY MODEL PROGRAMMER'S MANUAL 

The geometry module SIZE is the first operational subroutine to be called 
by MAIN, with NPC * 3. Input data passed from MAIN to SIZE by COMMON and NPC 
is changed to 2, and the subsequent call to SIZE results in geometric compu- 
tations which follow the sequence discussed in Part 1 of this volume. No 
other subroutines are called by SIZE apart from the "utility" program TPALT, 
whicn provides atmospheric data. Subroutine TFALT and its function have been 
described in Volume I. As shown in the subsequent flow chart, Figure I I. 3.1, 
there exist no iterative loops in SIZE, and for this reason the computational 
sequence is particularly simple. 

The equations employed in subroutine SIZE have been presented in Section 
II. 1. The computer symbols employed in the routine, together with a description 
of input and output from SIZE have been presented in Section II. 2. 
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FIGURE II. 3.1 SUBROUTINE SIZE 2 
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ELFFOELODN*HN+ELPC+LC 

ELBM-ELF-ELFFC 

SF«HC* (2 . 5* (ELODN*HN^ELPC)+3. 14*LC) 
SWBM>3. 14*ELBM»ELBM/BMLOD 
XARBM*SWBM/4. /BMLOD 


SW-WG/WGS 

B«SQRT(AR*SW) 

CRCLW*2.*SW/B/(1.+SLM) 

CBARW» (2 . *CRCLW/ 3 . ) * ( ( 1 . ♦SLM) - (SLM/ ( 1 . ♦SLM) ) ) 
TSWPLE* ( 1 . -SLM) / ( 1 . ♦SLM) /AR+TAN (DLMC4* . 017453) 
SWPLE*ATAN (TSWPLE) 

TSWPTE* 3 . * (SLM- 1 . ) / ( 1 . ♦SLM) /AR+TAN (DLMC4* . 017453) 
SWPTE*ATAN (TSWPTE) 

FHP*SWF/SQRT(2.) 

IIP*FHP*TSWPLE/2. 

CROOTW*CRCLW-HP+FHP*TSWPTE/2. 

VBARJ1*.85*(ELF*SWF**2./ (SW*CBARW))*CH1 
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~\J0_ 
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CBARIIT* (2. *CRCLirr/3. ) * ( ( 1 . *SL *H) - (SLMH/ (1 . ♦SLMH) ) ) 
V B AR V * . 3 36* ( E LF * HC * * 2 . / ( SW* B) ) *CV 1 
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4 VBARV* VBARVX j 
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FIGURE II. 1.1 SUBROUTINE SIZE 
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SVT»VBARV*SW*CV2 

BVT-SQRT(SVT*ARVT) 

CRCLVT»2. *SVT/BVT/ (1 . ♦SLMV) 

CBARVT- ( 2 . *CRC LVT/ 3 . )*((!. ♦SLMV) - (SLMV/ ( 1 . ♦SLMV) ) ) 

ELTH-VBARH*SW*CBARW/SHT 

ELTV“VBARV*SW*B/SVT 

DNAODBARN 

BD»S#F 

CMAOCBARW 

BL-ELF 

TRW1NG»SLM 


JCAEETPAI£(A^^ 


1 riHLP=PCAB-PAiMBC 1 


c 
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